Photon Added Detection 



CO 

o 
o 

CN 
>v 



O 

CO 



> 

in 

00 

in 
o 

CO 

o 

22 

Ok 
^— > ■ 

G ■ 

s : 



A.M. Brahczyk, 1 Tobias J. Osborne, 1 - 2 Alexei Gilchrist, and T.C. Ralph 1 

1 Centre for Quantum Computer Technology and Department of Physics, University of Queensland, QLD 4072, Australia 

^Department of Mathematics, University of Bristol, UK. 
(Dated: February 1, 2008) 

The production of conditional quantum states and quantum operations based on the result of 
measurement is now seen as a key tool in quantum information and metrology. We propose a new 
type of photon number detector. It functions non-deterministically, but when successful, it has high 
fidelity. The detector, which makes use of an n-photon auxiliary Fock state and high efficiency 
Homodyne detection, allows a tunable tradeoff between fidelity and probability. By sacrificing 
probability of operation, an excellent approximation to a photon number detector is achieved. 

PACS numbers: 42.50Dv 



I. INTRODUCTION 

In quantum theory, measurements encapsulate our ob- 
servation of nature. They are the link between the ab- 
stract machinery of the theory and its observational con- 
sequences. Because of this, it is not surprising that of- 
ten new measurement techniques and strategies can drive 
new applications. Moreover, the production of condi- 
tional quantum states and quantum operations based 
on the results of measurement is now seen as a key 
tool in realizing quantum information processing goals 
0,0- In optical schemes, conditional measurements pro- 
vide an effective nonlinearity that allows optical quantum 
gates to be fashioned Q, 0, 0, IE El , an d the creation of 
highly entangled states suitable for quantum metrology 

Often, however, the ideal measurements envisioned in 
theoretical proposals are not so easily realized experi- 
mentally. Linear optics quantum computation schemes 
such as in require high efficiency selective detectors 
(detectors able to distinguish between zero, one and sev- 
eral photons). The most promising detector candidate 
in this regard is the visible-light photon counter (VLPC) 
[ill IT^ | which has achieved efficiencies of the order of 
88%. Unfortunately these detectors require extreme op- 
erating conditions and suffer from high dark-count rates. 

In this manuscript we introduce the idea of a non- 
deterministic detector based on photon added detection 
(PAD), where we make use of high efficiency homo- 
dyne detection and mix the input state with an |n) 
Fock state prior to detection. This detector works non- 
deterministically, and there is an essential trade-off be- 
tween the probability that the detector works and the de- 
gree to which the detector functions as an n-Fock state 
projector. When the detector fails, this is clearly sig- 
nalled in the output. The essence of the detecting scheme 
is based on the observation that if we use homodyne 
detection and post-select within a narrow band of 2A 
around x = then the detection will only be sensitive 




FIG. 1: The probability density of getting a particular x value 
if we measure the X quadrature using homodyne detection. 
Results shown for various initial Fock states. 



to even photon numbers, see figure ^ By careful use of 
quantum interference, we can make the detector act like 
a projector onto a particular photon number. 

The structure of the paper is as follows. First we will 
introduce the scheme in general, then focus on the limit- 
ing case where A = to motivate its function. We then 
consider the effect of a finite A and discuss the trade 
off between probability of operation and fidelity. Finally, 
before concluding, we examine the effect of detector in- 
efficiencies in our scheme. 



II. THE SCHEME 

In order to characterise how well the detector functions 
we shall calculate the ability of the detector to pick out 
an appropriate state \a p ) from an entangled state of the 
form 
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when we measure mode b. The normalisation is A/"o = 
and the parameter w defines a window of states, 
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from which we want to pick out the central compo- 
nent. The reason for choosing this comparison is two- 
fold. Firstly we are interested in states precisely of the 
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FIG. 2: Quantum circuit describing our detector arrange- 
ment. 



above form where the states \a n ) represent multi-mode 
states which wc arc conditioning by detection and post- 
sclcction. Secondly, this approach provides an easily 
computable measure of how close to a \p)(p\ projector 
the detector functions in this context, since this approach 
reduces to a characterisation of state preparation |13| . 

With this characterisation in mind, consider the circuit 
in figure |21 We have some multi-mode state and we 
wish to condition the state of mode(s) a dependent on a 
photon number measurement on mode b. For simplicity 
consider only a single n-photon Fock state component in 
mode b, the general case is recovered through additivity, 
i.e. |V>) = A/o^ n \tp^). The input state is then some 
state \i/j( n >) — \a n ) a \n) b \p) c , where \a n ) is the associated 
component in mode a and mode c is initially in a p- 
photon Fock state. After interacting on a beam-splitter 
of reflectivity cos 2 (oj) and undergoing a phase shift A on 
mode b the output state is 
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where we have used the fact that the overlap between the 
quadrature amplitude eigenstates and the number states 
is given by 



H n {x e ) _ e -$ x *-i 



(4) 



and H n (x) is the Hermite polynomial of order n. We 
have chosen the convention that the 8 = quadrature 
operator can be written in terms of the mode operators 
as X = (a + at )/ y/2. Notice that the quadrature phase 
angles 9 and 4> are effectively not independent of A and 
that without loss of generality we can absorb those terms 
into A (so we will take A — 6 + — > A). For simplicity 
we shall also take (j) = and set the overall phase of 
this component to zero, and hence we can also drop the 
quadrature angle subscript on x and y. Now consider the 
case where we use a 50:50 beam-splitter so that u> = ir/4 
and we set A = tt/2. With these conditions equation (J3J) 
reduces to 



e -^(x 2 +y 2 ) e iirp 
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x H m+q (x)H n+p _ (m+q) (y) (6) 

To see how this detecting scheme is only sensitive to 
the p-Fock component we focus on the limiting case of 
A = next. 



where S and c 1 are the bosonic creation operators for 
modes b and c respectively, c = cos(cj), s = sin(cj), and 
finally we also have the usual binomial coefficients (") = 



(u — v) \v\ ' 

Modes b and c are now detected using separate bal- 
anced homodyne detectors. To an excellent approx- 
imation such detectors can be modeled as projectors 
onto small ranges of quadrature amplitude eigenstates 
\xg) where xg is a continuous variable with infinite di- 
mension, and 6 describes the phase relationship with 
the local oscillator of the homodyne detector. The fi- 
nal conditional state (unnormalised) , given we obtain 
xg in one detector and y^ in the other, is \ip C ond) — 



III. LIMITING CASE 

Consider only the special case where we happen to de- 
tect x = y = in the homodyne detectors. For these 
values, we can use 



n odd 

H n (0) = < 

I (n/2)I 



n even 



(7) 



This relation implies that only terms with even m + q 
will be non-zero, which in turn implies that n + p must 
be even also. If we now write g — > [g(n,p) + g'(n,p)]/2 
where g'(n,p) simply has the order of the summations 
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FIG. 3: Probability density plots for obtaining a particular x and y for the homodyne detections given an auxiliary photon 
number of |p) = |4), for various input number states \n). Only the x axis shown as the distributions are rotationally symmetric. 
By post-selecting on a narrow band near x = the detector becomes only sensitive to components with |n) = \p). Also, the 
noise form having a finite post-selection band comes from the nearby number states from the target state. 



reversed, we get 
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where we have set n = p+2k and used the fact that m+q 
must be even. From this expression it is clear that terms 
with odd k will also vanish. Terms with even k > will 
also vanish — this can be readily verified numerically. 
This then only leaves the terms with k = (n = p) as 
contributing to the state and so the detector picks 
out the \a p ) component. 

This analysis assumes an infinitesimal acceptance band 
for the detector. In order to assess the practicalities of 
the system we need to integrate over some range of val- 
ues around x = y = and evaluate success and failure 
probabilities. Clearly there will be a tradeoff between 
how well we project onto the p-photon Fock state and 
the probability of obtaining a successful outcome. 



IV. FINITE A 

The probability density for obtaining a value x in mode 
c and y in mode b will be 

P(x,y) = tr{|a;}<a;|<8>|y}<y|p} (9) 
= tr a {(x,y\p\x,y)} (10) 

where p is the three mode density matrix describing the 
state after the beam-splitter. This distribution is radi- 
ally symmetric about the origin, so we will switch to the 
polar co-ordinates r and 9 (where r 2 = x 2 + y 2 ) and ac- 
cept a particular result if it lies within a certain radius 
A. Intuitively we can see what the effect will be from 
figure 01 As we make A larger, the probability that a 
result falls within the accepted band, picks up contribu- 
tions from nearby states to the target state, and these 



will contribute to the error. The total probability that 
we get < r < A is 



P A = 2vr / P(r, 6)rdr 
/o 



(11) 



The (unnormalised) state immediately after destructively 
obtaining a particular x and y in the first two modes 
is Pa — (x,y\p\x,y). Consequently the ensemble of 
states that we would obtain if we where to only accept 
values within a radius A, would be 



Pa 



1 



d0 



' dr pW) 



(12) 
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FIG. 4: The main source of error for the detector is due to 
contributions from number states near the target state. Here 
we plot the difference in fidelity F, between two successive 
window sizes, w and w + 1. As can be seen, increasing the 
size of the window of states we are testing against makes little 
difference past a few states, consequently we will adopt w = 2 
in calculations. Note that A = 0.1 in the plots. 

To compare how well such a projector functions we can 
use the fidelity against the target state \a p ): 



F(A) = |(«>a|a p )| 



(13) 
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Note that in calculating this quantity we will assume that 



the 



are orthonormal. 



One of the important features of the PAD scheme is 
that it is sensitive only to a band of number states near 
the target state. This effect can be seen in the behaviour 
of the probability densities for states far away from the 
target state in figure [31 and is clearly demonstrated in 
figure where we show the rapid convergence in fidelity 
as we increase the number of nearby states to the one we 
are projecting out. 

As we increase A, the probability that we get a result 
we will accept also increases, but due to the overlap with 
the states near the target state the fidelity of the detector 
will drop. The actual probability is not a meaningful 
quantity in this context as it depends as much on the 
test state Q as on the parameters of the detector. The 
quantity we will use instead is a probability rate R = 
Pa/ '-Pidealj which is the probability we get divided by the 
expected probability if we had an ideal photo-counter. 
The tradeoff between fidelity and probability is quantified 
in figure [3] 




Visible-light photon counters can be modelled as ideal, 
but inefficient photon counters, at least for small photon 
numbers |T^ . 

The fidelity of the ideal detector in picking out the 
state \a p ) when used with the input state JTJ is then 



deal 



\(a p \Tr b {ILpPi n }\a p ) 




(15) 



where the summation extends to the maximum photon 
number, so for the test state in Q "max ~ p + w. 

For the PAD detector we can model inefficiencies sim- 
ply by considering a beam splitter of transitivity ry in 
front of both homodyne detectors 0] . The first observa- 
tion we make is that for high efficiency, the ideal detector 
obtains a higher fidelity. The trend with higher photon 
number is similar for both detectors. Where the advan- 
tage lies for the PAD is that the efficiency for current 
homodyne detectors is very high compared with avail- 
able photon counters. 

For a particular A and rj we can consider an equivalent 
ideal detector that gives the same fidelity. Constructing 
an equivalence in this fashion is particularly useful and 
was considered by where they compared an ideal pho- 
ton counter with homodyne detection in the context of 
quantum communication. As such, they used the mutual 
information as a means of comparison. For our scheme, 
we envision state preparation as the main application so 
we will use the fidelity as a means of comparison. This 
comparison is plotted in figure® for the ability to project 
out the state |ai) from the input state J2 n =o \ a n)\ n )- A 
detector able to achieve this projection forms a selective 
detector which is needed in many linear optics schemes. 



FIG. 5: The fidelity of operation for various target states \a p ) 
from the distributions Q is given (with w = 2). The curves 
are for fixed probability rates R. 



V. INEFFICIENT DETECTION 

The calculations so far have assumed unit efficiency de- 
tection. In this section we explore the effect of non-unit 
detection efficiencies for the PAD, although it should be 
noted from the outset that detection efficiency for homo- 
dyne detection is very high (in the region of 98% 0). 
We will compare the performance of the PAD to an ideal, 
but inefficient photon counter, which we model by the 
POVM elements n p : p = 0, 1, . . ., where p is the number 
of detected photons, with 



ni—p 



rf(l-n) m -*\m)(m\ 



(14) 




FIG. 6: Equivalent ideal single photon detector efficiency as 
a function of the acceptance width A, and the Homodyne 
efficiency rj. 
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VI. DISCUSSION AND CONCLUSIONS 

Because of it's non-deterministic nature, we envision 
applications of this detector mainly in state preparation, 
where non-classical states are prepared through condi- 
tioning on photon number detection. We could prepare 
a good approximation to an \n) photon state required by 
our detector, by using spontaneous parametric down con- 
version and a detector cascade in one arm. Even if the 
detectors in the cascade are inefficient, if, say three detec- 
tors register a click, then we have at least a three photon 
term in the other arm. The errors caused by having more 
than the required number of photons are offset by the low 
probability of such events. One intriguing possibility is 
to employ this detector in a proposal by Dakna, et at 
|l8j |. In the Dakna scheme, a good approximation to an 
optical Schrodinder cat state is generated by mixing a 
single mode squeezed state on a beam-splitter with the 
vacuum and conditioning on detecting a certain number 



of photons in one of the exit ports. 

Another possible extension is to use other parameter 
choices, and post-selection choices to directly project out 
certain distributions of photon number terms. 

We have presented a non-deterministic scheme which 
functions as a high-fidelity Fock state projector. This de- 
tecting scheme allows a tunable tradeoff between the fi- 
delity and probability of detection. The weaknesses of the 
scheme are that it requires an |n) photon state and that 
it is non-deterministic. The |n) photon state could be 
prepared in the first instance simply by conditioning the 
output of a spontaneous parametric down converter with 
a traditional detector cascade. The non-deterministic na- 
ture of the scheme leads us to conclude that the main 
application for the detector will be in state generation. 
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tralian Research Council. 
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